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Dipolar interactions support the formation of inter-site soliton molecules in a stack of quasi- 
1D traps. We show that the stability and properties of individual solitons, and soliton molecules 
in such a geometry crucially depend on the interplay between contact and dipolar interactions. 
In particular, two different quasi- ID soliton regimes are possible: a ID soliton characterized by 
purely repulsive DDI and a 3D soliton for which a sufficiently large dipole moment renders the DDI 
attractive. Furthermore, we find that contrary to the case of dimers of polar molecules, the soliton 
dimers exhibit a nontrivial behavior of the elementary excitations that stems from the competition 
between on-site and inter-site DDI. Finally, we prove the existence of soliton trimers in a regime 
where molecular trimers do not occur. We demonstrate that the soliton molecules that we report 
are well feasible under realistic experimental conditions. 

PACS numbers: 03.75.Lm, 03.75.Kk, 05.30.Jp 



I. INTRODUCTION 

Recent developments in experiments on ultra-cold po- 
lar molecules [1-3], atoms with large magnetic dipole mo- 
ment [4-7] , and Rydberg atoms [8] opens new promising 
perspectives in the rapidly progressing research on dipo- 
lar quantum gases. Interestingly, the presence of long- 
range and anisotropic dipole-dipole interactions (DDI) 
essentially modifies the behavior of quantum gases lead- 
ing to a wealth of a new physics [9, 10]. 

Dipolar effects are particularly relevant to what con- 
cerns the nonlinear properties of dipolar Bose-Einstein 
condensates (BECs). Crucially, whereas non-dipolar 
BECs present a local Kerr-like type of nonlinearity, the 
nonlinearity in dipolar BECs exhibits a nonlocal charac- 
ter, similar to that in plasmas [11], photorefractive media 
[12, 13] and nematic liquid crystals [14, 15]. Interestingly, 
this nonlocality results in novel physical phenomena, in- 
cluding stabilization of two-dimensional solitons [16, 17]. 

The long-range character of the DDI plays a substan- 
tial role in the physics of dipolar gases in optical lattices, 
even in the absence of inter-site hopping. While a non- 
dipolar gas in such a deep lattice may be considered as a 
system of mutually independent gases, nonlocal inter-site 
interactions in a dipolar gas couple the disjoint sites. In 
particular, in the physics of polar molecules this feature 
gives rise to a variety of unprecedented few-body bound 
states such as inter-site dimers [18, 19], trimers [20, 21] 
and filaments [22, 23]. 

The inter-site interactions play also a key role in the 
behavior of a dipolar condensate in an optical lattice. 
Specifically, they have been found to fundamentally mod- 
ify the BEC excitation spectrum [24, 25] and to affect sig- 
nificantly the stability of the condensate [26, 27]. More- 
over, inter-site interactions may lead to a correlated mod- 
ulational instability, in which a locked density modula- 



tion pattern is shared among non-overlapping sites, after 
a quench of a condensate into instability. Interestingly, 
such correlated modulational instability may result in the 
dynamical formation of soliton filaments and crystalline 
structures [28]. 

In this paper we analyze in detail the physics of dipolar 
bright solitons in a stack of quasi- ID condensates. We 
focus on the stability and properties of soliton dimers 
and trimers, which constitute the building blocks of the 
above-mentioned soliton filaments and crystals, respec- 
tively. These two- and three-soliton bound states are an 
example of the so-called soliton molecules. Recently, an 
optical equivalent of such objects has been realized ex- 
perimentally in optical fibers [29, 30] and a variety of 
theoretical proposals to create atomic soliton molecules 
have been presented [31-33]. Soliton dimers share some 
properties with molecular dimers. However, as we dis- 
cuss in detail below, intra-soliton interactions (of course 
absent in the case of individual polar molecules) are deci- 
sive for their stability and elementary excitations. More- 
over, whereas molecular trimers may be found (in absence 
of any additional lattice [21]) only for a rather narrow 
window of the dipole moment orientations [20], soliton 
trimers may exist for the orientations for which trimers 
of individual polar molecules are precluded. 

The article is structured as follows. Sec. II introduces 
the general formalism. In Sec. Ill we compute the uni- 
versal stability diagram for a single dipolar soliton in a 
quasi- ID trap and we show that such geometry supports 
two stable soliton regimes differing substantially in the 
character of the dipolar interactions. Section IV is de- 
voted to the study of properties of the soliton dimers. 
We discuss the inter-soliton binding potential and the 
nontrivial dependence of the dimer elementary excita- 
tions on the dipolar coupling. In Sec. V we analyze the 
trimer case, showing that soliton trimers may be found 
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FIG. 1. (Color online) Scheme of a stack of quasi- ID tubes 
of dipolar Bose-Einstein condensates. 



in a regime where molecular trimers would be unstable. 
We conclude in section VI. 



II. MODEL 

In the following we consider a dipolar BEC loaded in a 
stack of M parallel quasi- ID traps (tubes), formed by a 
2D optical lattice with sites located at i/j = jA (Fig. 1). 
The inter-tube potential barrier is considered sufficiently 
large to suppress any hopping. In each tube we as- 
sume a strong harmonic confinement of frequency cj ± in 
the xy plane and no confinement along the z direction. 
The atoms possess a magnetic dipole moment \i (the re- 
sults are equally valid for electric dipoles, such as polar 
molecules) oriented along the y axis, in the side-by-side 
configuration, by a sufficiently large external field. Intro- 
ducing a wavefunction ^ (r) that describes an atomic 
cloud in a site j holding N atoms, the system of non- 
local coupled Gross-Pitaevskii equations (GPE) reads 
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Here, Uj (r) 

QdN (l — 3 cos 2 6) j |r — r'| 3 is the dipole-dipole poten- 
tial where g d = hqh 2 /An with //q being the vacuum per- 
meability and 6 the angle between the vectors joining 
two interacting particles and the direction of the dipole 
moment. The short-range interactions are characterized 
by g = 4:7ra sc h 2 /m with a sc being the s-wave scattering 
length. In the following we consider attractive short- 
range interactions (a sc < 0). 



III. DIPOLAR SOLITON IN A SINGLE 
QUASI-1D TRAP 

We discuss first the conditions of existence of a stable 
bright soliton in a single quasi- ID trap (M = 1). To this 
end we assume a 3D anisotropic Gaussian ansatz 
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where l x , l y , and l z are the variational widths along x, 
y and z directions, respectively. Employing this ansatz 
into Eq. (1) we obtain the energy of the system 
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that in the cases of our interest may be evaluated analyt- 
ically in terms of elliptic integrals [34]. A stable soliton 
solution corresponds to a minimum in the energy func- 
tional E (l x ,l y ,l z ) at finite non-zero values of the soli- 
ton widths. In Fig. 2 we present the universal stability 
diagram as a function of the dimensionless parameters 
= gNfahuJl, and g* = g d N/27rhujJ 3 ± . 
Interestingly, two different soliton regimes may be 
found, which differ remarkably in their properties and 
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FIG. 2. (Color online) Universal stability diagram for a dipo- 
lar bright soliton in a single quasi- ID trap. Three regimes 
occur: stable soliton, instability against a 3D collapse and 
soliton expansion along the axis of the trap. The dashed 
line represents the stability boundary for a soliton dimer with 
A = 6l ± . 
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FIG. 3. (Color online) Reentrant character of the soliton sta- 
bility in a single quasi- ID trap in the vicinity of g*. Here, 
|<f I = 0.95. 



stability for growing g d > 0. For sufficiently small 
l#*l < \dc\i with \g*\ ~ 1, a soliton may be considered 
as purely ID, i.e. l x 



l y ~ l ± = y^h/muj i , whereas 



l z ^> l±. For such soliton, the DDI remains repulsive for 
any g d . As a result, the soliton width l z increases mono- 
tonically for growing g^ , until diverging at a critical value 
at which the soliton delocalizes. The condition for soliton 
stability against the expansion may be then found analyt- 
ically from Eq. (3), \g*\/g* d > 2tt/3 (straight solid line in 
Fig. 1). On the contrary, for \g*\ > \g*\ the atomic cloud 
cannot be considered any more as ID, since l z becomes 
comparable with the transversal widths. As a result, a 
stable soliton solution occurs that clearly displays a 3D 
character. In this regime, the DDI interaction changes its 
character from repulsive to attractive at a finite g* d > 
value, and hence for further growing g* d the soliton width 
decreases until the soliton becomes unstable against 3D 
collapse. Furthermore, we note that in the vicinity of 
|g*|, the stability diagram presents an interesting reen- 
trant character as a function of g dl first expanding, then 
re-binding and finally collapsing (Fig. 3). Interestingly, 
contrary to the soliton-expansion transition, at which the 
soliton width smoothly diverges, the re-binding transi- 
tion is first-order-like, since the soliton abruptly re-binds 
at a finite width. 



IV. SOLITON DIMERS 

We assume in the following that a soliton in each tube 
is in the ID regime discussed in Sec. Ill (this condition 
is self-consistently verified). At the end of this section 
we briefly comment on the case of solitons in the 3D 
regime. In the ID regime, the wavefunctions factor ize 
\&j (r) = <j)j (x,y) ipj (z), with <pj (x,y) the ground state 
wave function of the transverse harmonic oscillator in 
a site j. Employing the convolution theorem [35] and 
integrating Eq. (1) with respect to x and y, we arrive at 



the dimensionally reduced system of equations 
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with n m (k z ) the Fourier transform of the axial wave func- 
tion density n m (z) = \ip m (z)\ 2 in a site m and 
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For stable individual solitons the inter-site DDI may re- 
sult for g d > in a binding of two solitons in different 
quasi- ID tubes into a soliton dimer (Fig. 4). This dimer 
resembles the case of the recently reported dimers of in- 
dividual polar molecules. However, as discussed below, 
the interplay between intra-soliton interactions and inter- 
soliton interactions leads to a non-trivial effects in the 
physics of the soliton dimer, which do not occur in the 
case of molecular dimers due to the absence of on-site 
DDI. 

Two solitons localized in neighboring quasi- ID tubes 
(j = 0, 1) and with a relative displacement z r along the 
axis direction z (Fig. 4), experience an interaction poten- 
tial 
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FIG. 4. (Color online) Inter-soliton binding potential for the 
case of the soliton dimer. The red dashed line represents 
the potential calculated within the point-like approximation 
E%. The blue solid line shows the actual potential com- 
puted numerically with Eq. (7). Here, we consider the case 
of 52 Cr condensate (/x = 6/xs, with \ib the Bohr magneton), 
dsc = — 7.1 ao (with do the Bohr radius), N = 100, A = 
6/^ = 512 nm, and the lattice potential depth s = 13.3Er 
(recoil energy). These parameters refer to uj±_ — 26.7 kHz 
and (<7*,<7d) = (—0.88,0.45). The inset depicts schematically 
the soliton dimer arrangement. 



We calculate E]j(z r ) evolving Eq. (5) in imaginary time 
to obtain the ground state of the dimer ipj(z) and 
then shifting the solitons to the distance z r . Due to 
the anisotropy of the DDI the inter-soliton potential is 
maximally attractive for z r = 0, becoming repulsive 
for large z r (Fig. 4). Naturally, the binding potential 
E]j(z r ) calculated for actual soliton wave-packets is sig- 
nificantly weaker than that expected for point-like parti- 
cles E° D (z r ) = g d N(z 2 - 2A 2 )/(z 2 + A 2 ) 5 / 2 . Neverthe- 
less, we note that even for the case of the relatively small 
dipole moment of 52 Cr, which we employed in our calcu- 
lations for Fig. 4, the energy scale of the binding remains 
significant (~ 100 Hz). The binding would be of course 
stronger for condensates of atoms with larger dipole mo- 
ment, such as dysprosium [6] and erbium [7], or in the 
case of polar molecules [1-3]. 

We now focus on the essential properties of the soliton 
dimer. First, following the imaginary time evolution of 
Eq. (5), for a given A/Zj_, we compute the width l z of the 
solitons forming the dimer as a function of g* and (see 
Fig. 5 (top)). Since we consider the ID soliton regime, 
with an overall repulsive intra-soliton DDI, an increase of 
gd results in a broadening of the solitons, and eventually 
to the instability of the individual solitons against ex- 
pansion. Note, however, that the attractive inter-soliton 
interactions, while providing the binding mechanism it- 
self, induces a trapping of each soliton by its neighbor, 
which contributes to stabilization of each soliton against 
expansion. This increases the stability threshold found in 
Sec. Ill for an individual soliton, as shown by the straight 
dashed line (at \g*\/g^ — 1-78) in Fig. 2, obtained from a 
similar 3D variational calculation as that of the previous 
section. 

The properties of the soliton dimer must be compared 
with those of inter-site dimers formed by individual po- 
lar molecules. In the latter case, the localization of 
each molecular wave-packet is solely due to the attractive 
inter-site DDI, which induce a mutual trapping of both 
molecules. This means, in particular, that for g* d = 
each of the wave-packets delocalizes. Furthermore, ow- 
ing to the absence of intra-wave-packet repulsive DDI, 
an increase of g d can only amplify localization and so 
the molecular dimer width decreases monotonically as a 
function of g^, unlike the case of the soliton dimer. As 
a result, molecular dimers become stiffer (i.e. present 
growing excitation energies) for growing DDI. 

On the contrary, the lowest-lying excitation of the soli- 
ton dimer presents a more involved behavior due to the 
interplay between intra- and inter-soliton DDI. We study 
the lowest-lying excitations by monitoring the real-time 
dynamics of the solitons following a small distortion of 
the ground state solution in the form ipj(x,t = 0) = 



(0) -{(kjx+PjX 2 ) 



, corresponding to a perturbation of the 



soliton positions and their widths. Fig. 5 (bottom) shows 
the result of the Fourier transform of the position (z(t)) 
of one of the two oscillating solitons and hence the fre- 
quency of the dimer lowest-lying excitation (this is veri- 
fied additionally by inspecting the Fourier transforms of 
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FIG. 5. (Color online) (top) Width of the soliton dimer as 
a function of g* d for \g*\ = 0.6 (blue dashed line) and \g*\ = 
0.8 (green dot-dashed line), for the same parameters as in 
Fig. 4. The vertical dashed lines indicate the dimer expansion 
threshold, (bottom) Frequency of elementary excitations of 
the soliton dimer for the same parameters. The inset shows 
a scheme of the dimer elementary excitation mode. 



soliton width and density oscillations). For sufficiently 
small DDI, and so for a small solitons widths, the lowest- 
lying excited mode of the dimer is associated exclusively 
to the motion of the center-of-mass of each soliton. In 
consequence, as g d grows, so does the energy of dimer 
excitations, resembling the case of molecular dimers. In 
contrast to the molecular dimers, however, after reaching 
a certain critical value of g d the soliton dimer becomes 
progressively softer (i.e. it exhibits decreasing excitation 
energies). This phenomenon arises because the soliton 
widths increase due to the repulsive intra-soliton DDI, 
and as a result the lowest-lying excitation becomes even- 
tually an admixture of both position and width distor- 
tions. As discussed before, for a sufficiently large g d the 
dimer becomes eventually unstable against expansion. 

Finally, we stress that soliton dimers may exist as well 
in the 3D regime defined in Sec. Ill, i.e. for |g*| > \g*\. 
As depicted in Fig. 2, the stability threshold against the 
soliton dimer collapse is basically the same as that for an 
individual soliton. Contrary to the ID case, in the 3D 
regime the width of a soliton is relatively small l z ~ l ± 
and so the binding potential between the two solitons, 
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such as the one depicted in Fig. 4, becomes comparably 
deeper, approaching the point-like approximation E® D . 
Moreover, for |g*| > |g*| the soliton width never becomes 
large enough to cause the mixing of position and width 
excitations. As a result, in the 3D regime, for growing 
gd values the soliton dimer becomes only stiffer, up to 
the collapse threshold, similary to the case of molecular 
dimer s. 



V. SOLITON TRIMERS 

Interestingly, the DDI may lead to the formation of 
soliton molecules comprising of more than two solitons, in 
particular soliton trimers (Fig. 6). We note that trimers 
(and even more involved complexes) have been predicted 
as well for individual polar molecules [20, 21]. How- 
ever, molecular trimers have been found to exist only 
in a rather narrow window of dipole moment orienta- 
tions with respect to the trap axis, in the very vicinity 
of the magic angle, such that intra-site repulsion is min- 
imized and inter-site attraction is maximized. In partic- 
ular, molecular trimers are precluded if the dipole ori- 
entation is aligned along the trap axis. Furthermore, as 
noted in Sec. IV, the formation of molecular bound states 
is handicapped by the fact that the inter-site interactions 
do not only provide a binding between the molecules but 
are also indispensable for the localization of the individ- 
ual molecular wave-packets themselves. This contrasts 
with the soliton case, where the existence of localized 
wavepackets is supported by intra-soliton interactions. 
As a result, as we discuss in this section, the interplay 
between inter- and intra-soliton interactions allows for 
stable soliton trimers for dipole moment orientations in 
which molecular trimers are absent. 
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FIG. 6. Scheme of the soliton trimer. The dipole moments 
are aligned in the head-to-tail configuration providing attrac- 
tive intra-site and repulsive inter-site dipolar interactions. In 
our work we mimick this scheme with qualitatively equivalent 
arrangement of dipoles aligned along the y axis (side-by-side 
configuration) but with < (see text). 
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FIG. 7. (Color online) Universal stability diagram of a single 
dipolar bright soliton in a quasi- ID trap for g^ < 0. Two 
regimes occur: stable soliton and 3D instability against col- 
lapse. We indicate additionally the regime of the trimer in- 
stability against fusion (see discussion in text). 



In the following we consider for theoretical simplic- 
ity the case of dipoles oriented along the y axis (in the 
side-by-side configuration as that of the soliton dimer) 
but with gd < 0. This may be achieved by means 
of a rotating magnetic field [36], or microwave dress- 
ing for polar molecules [37]. The results would be how- 
ever qualitatively very similar to the case of dipoles ori- 
ented along the tubes, since both cases are character- 
ized by repulsive inter-site DDI and attractive intra-site 
DDI. Although the attractive inter-site interactions seem 
naively to involve soliton fusion in the bottom tube, and 
hence to preclude the existence of the soliton trimer, 
such trimer results actually from a nontrivial interplay 
between inter-tube repulsion and intra-tube attraction. 
Namely, the single soliton in the upper tube provides 
a repulsive potential barrier that prevents the fusion of 
the two mutually-attracting solitons in the bottom tube, 
hence keeping the soliton trimer stable. 

A major difference with respect to soliton dimers lies 
in the fact that now g^ < 0, and hence the intra-soliton 
interaction is attractive. As a result for growing \g^\ the 
individual solitons shrink, i.e. the trimer is not unsta- 
ble against the expansion of the individual solitons but 
rather against their collapse, since the solitons become 
eventually 3D for a sufficiently large \g^\. As shown in 
the dimer case, the threshold for the collapse instabil- 
ity is basically given by the intra-soliton physics. We 
have hence analyzed the stability of a soliton in a single 
quasi- ID trap for gd < (see Fig. 7), using the same 3D 
variational Gaussian ansatz discussed in Sec. III. Natu- 
rally, soliton trimers may exist only within the stability 
region of individual solitons. 

In the following we analyze the properties of trimers 
well within the ID regime, i.e. far from the 3D collapse 
threshold, for which we can safely employ the ID GPEs 
(Eq. (5)). In particular, after obtaining the ground state 
of the trimer configuration by means of the imaginary 
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time evolution of these equations, we have computed the 
binding potential of the trimer Et(t) (Fig. 8) as a func- 
tion of the distance r between the solitons in the bottom 
tube (Fig. 6). Crucially, at an intermediate distance r m i n 
Et{t) shows a local minimum that offers the possibility 
of a soliton trimer. A point-like approximation of the 
solitons would induce a binding 



E° T (r) = g d N 



1 16(r 2 - 8A 2 ) 



(r 2 + 4A 2 ) 5 / 2 



(8) 



resulting in an equilibrium position r^ in /A ~ 3.73, inde- 
pendently of ga- This approximation, however, departs 
significantly from the actual binding potential Er(r), 
proving again the relevance of the spatial extension of 
solitons. Specifically, as shown in Fig. 9 (top), the trimer 
size, understood as the actual equilibrium distance r min , 
decreases with growing \g^\ (whereas the binding energy 
increases). We also note that, as it may be expected, 
the soliton trimer is more loosely bound than the soliton 
dimer. For typical parameters of 52 Cr condensate, which 
we employed for the Fig. 8, the binding energy is of the 
order of 10Hz. We stress, however, that the binding will 
be certainly stronger in the case of more magnetic atoms 
(Dy, Er) or polar molecules. 

Note that contrary to the dimer case, the soliton trimer 
is related with a local minimum of the energy functional. 
In particular, the global energy minimum results from 
the fusion of the solitons in the bottom tube into a sin- 
gle soliton, which forms a tilted dimer with the top soli- 
ton. The trimer configuration of Fig. 6 is hence strictly 
speaking a metastable solution, which is separated from 
the fused solution by a potential barrier (Fig. 8). Macro- 
scopic quantum tunneling through this barrier is negligi- 
ble, and hence the metastable solution may be considered 




FIG. 8. (Color online) Soliton trimer potential energy. The 
red dashed line represents the potential obtained from 
the point-like approximation, whereas the blue solid line de- 
picts the actual binding potential Et calculated numerically 
integrating Eqs. (5). We consider 52 Cr BEC (with g d < 0) 
with a sc — — 4.0 ao and N — 100 atoms in every soliton, i.e. 
(9*,9d) = (-0.50,-0.45). Here, A = 512 nm, s = 13.3E R 
and uj±_ — 26.7 kHz. The inset shows the potential energy 
minimum which sustains the trimer bound state. 
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FIG. 9. (Color online) (top) Size of the soliton trimer as a 
function of g* d for \g*\ = 0.25 and the remaining lattice param- 
eters as that in the Fig. 8. The dashed line indicates the result 
from the point-like approximation E^. (bottom) Trimer ele- 
mentary excitations for the same parameters. The insets show 
schemes of the trimer elementary excitation modes. 



for all practical purposes as stable (as we have checked 
in real-time evolution). The potential barrier dissappears 
at a sufficiently small \g^\, at which the soliton trimer be- 
comes abruptly unstable against soliton fusion (see Fig. 7 
and Fig. 9). 

Finally, as for the case of the soliton dimer s, we 
have analyzed the lowest-lying excitations of the soliton 
trimer. Since now gd < 0, the solitons are always well 
localized. Hence, contrary to the dimer case, the lowest 
lying excitations are related solely to the solitons center 
of mass motion (without an excitation of the width of 
the solitons). We may hence define two different types of 
elementary excitations, characterized by an in-phase and 
an out-of-phase motion of the soliton pair in the bottom 
tube, respectively (Fig. 9 (bottom)). As in Sec. IV, we 
have probed these modes perturbing the soliton widths 
and positions of the trimer ground state, and monitoring 
the subsequent real-time dynamics given by Eq. (5). Af- 
ter Fourier transforming the soliton positions, we obtain 
the lowest-lying excitations as a function of \g^\. The 
results of the two excitation frequencies are depicted in 
Fig. 9 (bottom), which shows that for all \g^\ the out- 
of-phase mode is always less energetic than the in-phase 
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mode. 



VI. CONCLUSIONS 

In summary, inter-site dipolar interactions support the 
formation of soliton molecules in a stack of quasi- ID 
tubes. The stability properties of quasi- ID solitons and 
inter-site soliton molecules depend crucially on the in- 
terplay between dipolar and contact interactions, and 
the competition between intra-site and inter-site effects. 
In particular, two different quasi- ID soliton regimes are 
possible: ID solitons for which the intra-soliton DDI is 
always repulsive and that become eventually unstable 
against soliton derealization, and 3D solitons, for which 
the DDI changes its character from repulsive to attractive 
for growing DDI, and that become eventually unstable 
against soliton collapse. We have shown that, contrary to 
the case of dimers of individual polar molecules, the inter- 
play between intra-soliton interactions and inter-soliton 
DDI leads to a non-trivial behavior of the lowest-lying 
excitations of soliton dimers. In the purely- ID regime 
a growing DDI render the dimer stiffer up to a maxi- 



mum beyond which an increasing DDI softens the dimer 
due to the admixture between position and width excita- 
tions. Finally, we have shown that soliton trimers may be 
constructed for attractive intra-site and repulsive inter- 
site DDI due to a subtle interplay between intra-tube 
attraction and inter-tube repulsion. Interestingly, these 
trimers occur in a regime in which trimers of individ- 
ual polar molecules are not possible. The reported soli- 
ton molecules can be observed under realistic conditions 
within current experimental feasibilities. Moreover, we 
emphasize that the soliton binding mechanism described 
in this work can be straightforwardly generalized to en- 
gineer even more intricate soliton complexes comprising 
a larger number of solitons in more sites of an optical 
lattice. 
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